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Abstract 



On 3 and 4 October, 1995, we obtained high spectral resolving power (R= A/AA = 60,000) 
observations of comet 122P/de Vico using the 2DCoude cross-dispersed echelle spectrograph on the 
2.7-m telescope of McDonald Observatory. The spectra cover the wavelength range from 3830A- 
10192A. The spectra fr om 3830-5776A are continuous; from 5777-10192A there are increasing 
interorder gaps. The comet was at a heliocentric distance of 0.66 AU and a geocentric distance of 
1.0 AU. Comet de Vico has a very high gas-to-dust ratio and the spectra have excellent signal/noise. 
These two factors combined to yield spectra with a large number of emission lines. We have collected 
laboratory molecular line lists and have used these line lists in order to identify as many of the 
detected lines as possible. We have identified 12,219 emission lines and have located another 4,055 
lines which we cannot identify. We present representative spectra and identifications along with a 
description of our plans to make these spectra and identifications available to the community. This 
atlas should prove a valuable tool for future studies of comets. 
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1 Introduction 

Understanding the physics and chemistry of comets can lead to significant improvements in our 
knowledge of the formation of the Solar System. Such studies allow constraints on the physical 
conditions in the solar nebula, on the dynamical evolution of the disk and of planetesimals, and 
on the chemical evolution of these bodies. High resolving power spectroscopic observations offer 
unique opportunities for studies of the origin and subsequent chemistry of comets. For example, 
measurement of atomic isotope ratios in a sample of "Kuiper belt" and "Oort cloud" comets can 
place valuable new constraints on physical processes in the solar nebula and proto-planetary disk 
during the comet formation epoch. 

The complex molecular spectrum of the coma of a comet results in thousands of lines being 
visible in the spectral region from 3800-10200A. In order to undertake studies of a particular 
feature or molecule, one must be able to identify most of the lines in the region of interest. Yet 
many spectral regions are complicated by the superposition of lines from more than one molecule. 
In addition, laboratory studies of diatomic and polyatomic molecules are often difficult to obtain 
and incomplete. Valk et al. (1992) published a catalogue of line identifications for the near-UV 
and blue region of the spectrum of Comet C/1989X1 (Austin) from spectra at moderate spectral 
resolving power (R=2,500). Brown et al. (1996) published a catalog of cometary emission lines 
from comets 109P/Swift-Tuttle and C/1989X1 (Austin). They found 2997 lines and were able to 
identify 2438 of them. Zhang et al. (2001) utilized high resolution spectra of C/1995 01 (Hale- 
Bopp) and the catalogue of Brown et al. in order to identify 532 emission lines (73 unidentified). 
These catalogues are a critical first step for cometary science. 

This paper represents an update to the catalogues of Brown et al. and Zhang et al. We 
present observations of comet 122P/de Vico obtained at McDonald Observatory in October 1995. 
Comet de Vico is the perfect target for such a study since it has a very high gas-to-dust ratio 
and was reasonably bright, allowing for spectra of exceptionally high signal/noise with very little 
continuum (signal/noise at peak was >300). We identified 12,219 lines in the spectrum of de Vico. 
In addition, there are 4,055 unidentified lines in our spectra. 

2 Observations and Reductions 

We observed comet 122P/de Vico with the 2.7-m Harlan Smith telescope of McDonald Observatory 
and the 2DCoude spectrograph (Tull et al. 1995) during October 1995 (see table | for information 
on the observations). These observations were obtained at a spectral resolving power of 60,000 
and covered the wavelength range from 3800-5776A continuously, with spectra with increasing 
interorder gaps from 5776-10192A. Table [IJ indicates the spectral coverage of each of the echelle 
orders. In all cases, the slit was 1.2 arcsec wide by 8.2 arcsec long. The spectra were extracted 
using variance weighting along the spatial dimension, with no attempt made to preserve spatial 
information. The signal/noise was extremely high throughout. De Vico represents a perfect comet 
for such identification work since the gas/dust ratio is one of the highest known for a comet. 
This results in spectra with almost no continuum so that even weak features are unambiguously 
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detected. In addition to spectra on the optocenter, we obtained one spectrum with the telescope 
offset 100 arcsec into the tail. These tail spectra are useful for distinguishing between neutral and 
ionized species. 

Incandescent lamp spectra were obtained in order to flat field the data. The wavelength 
calibration was made by fitting a two-dimensional function to observations of a ThAr hollow- 
cathode lamp. Over 3100 lines were fit, with a fourth order fit along an order and fifth order 
fit across orders, resulting in an rms error of the wavelength of ~2.5mA. This is comparable to, 
or better than, typical laboratory spectra described in this paper. The cometary spectra were 
Doppler-shifted by a velocity equal to the geocentric radial velocity of the comet (i.e by —A from 
table | to bring them onto the laboratory rest frame). 

A cometary spectrum consists of molecular and atomic emission lines, superimposed on 
a continuum which results from sunlight reflecting off the dust. It is necessary to remove the 
underlying continuum from the emission spectrum if the goal of a study of the spectra is to determine 
the amount of a particular gas. However, removal of the continuum is fraught with uncertainties in 
the color of the dust and in photon statistics. Thus, removal of the continuum can induce noise into 
spectra. Fortunately, comet de Vico has a high gas/dust and the continuum is essentially negligible 
at all wavelengths. Given this and the problems associated with continuum removal, we chose not 
to remove the continuum for our purposes in this paper of line identification. This results in an 
offset of the spectrum from zero (typically the continuum for de Vico was about 350 counts and 
the emission lines are many thousands of counts). The continuum can also result in the apparent 
diminution of the strengths of emission lines when they coincide with a solar absorption feature. 
This might result in our not detecting a very weak feature. 

In addition to the underlying continuum spectrum, there is a superposition of a telluric 
spectrum onto the cometary spectrum. This spectrum consists of emissions due to metastable 
oxygen, sodium, OH and O2 and absorption due to O2 and H2O. The telluric spectrum is at the 
rest frame of the Earth's atmosphere and is, therefore, shifted by —A with respect to the cometary 
rest frame. Thus, we Doppler shifted by this same velocity the telluric emission line lists discussed 
in the next section to identify the telluric features present in the Doppler-shifted cometary spectra. 

The telluric emission lines are generally narrow and only affect the wavelength at which 
they occur. The telluric absorption is a different problem entirely. In the regions of the strong 
absorptions, such as the O2 "A" band or the 9300A H2O band, almost all photons entering the 
atmosphere are absorbed before they reach the telescope. Thus, even if a line at these wavelengths 
is excited in the cometary spectrum, we will not record the feature in our spectrograph. There is 
no way to remove the effect of these saturated telluric absorptions and the information from the 
comet is forever lost (unless a spectrum taken at a different time has a Doppler shift that moves 
the cometary line and the telluric absorption apart). The orders most affected by the telluric 
absorption are orders 50 (0 2 "B" band), 48 and 47 (H 2 "a" band), 45 (0 2 "A" band), 42 and 41 
(H 2 "z" band), and 38-35 (H 2 "7" and "p" bands). 
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3 Linelists 

One of the most difficult aspects of compiling a cometary spectral atlas is finding good laboratory 
line lists of all of the relevant molecules. As a first step in building our spectral atlas, we searched 
out many molecular line lists. Many of the molecules of importance were studied years ago and 
the line lists do not exist in digital form; we have spent many hours scanning line lists and typing 
in and checking large quantities of numbers by hand (possibly with undetected errors introduced). 
Some of the relevant molecules are not well understood and there are many unassigned lines. For 
some molecules, we have used work by several different authors to be as complete as possible. In 
the case of some molecules, such as CN, we had access to line lists which were much more inclusive 
than necessary. The tables of the identified molecular lines which we include as part of this paper 
list only those lines which we believe we detected. However, such over-complete line lists may prove 
useful for future cometary spectra. We note, also, that in the red and near-IR, there are interorder 
gaps in our data. In these regions, we leave gaps in the lists of identified lines. Refer to table |TJ for 
a complete list of the wavelengths for which we have observations. 

3.1 C 2 

There are two principal band systems of C2 which are observed in the optical in the spectra of 
comets. These are the Swan, or d 3 II 9 - a 3 ILj, system, and the Phillips, or A 1 H U - X 1 ^^", system. 
The Swan system is dominant in the green, orange and red region of the spectrum, while the 
Phillips is important in the near-IR and IR. 

For the Swan system, we used the line list of Phillips and Davis (1968). We extended this 
line list for high J-values of the 0-0, 0-1 and 2-1 bands using the wavelength spacing of the highest 
previously defined lines to identify which emission lines were the continuation of these bands. We 
derived slightly different wavelengths for four of the highest J- value lines of the R-branches of the 1- 
band than in the orig inal reference (AA =0.042-0.0761). Most of our extensions are in agreement 
with the additional wavelengths also derived by Brown et al. (1996), though we observed to higher 
J-levels for the 0-0 band. 

For the Phillips system, we used the works of Ballik and Ramsay (1967) and Chauville et al. 
(1977). 

3.2 C 3 

The C3 molecule was first observed in cometary spectra in 1882, although the species causing the 
emission feature in comets was unknown. It is difficult to study in the laboratory because it is 
relatively unstable. Thus, the molecule causing the emission lines was not identified until 1951. In 
cometary spectra, the important bands are known as the "4050-A Group" , indicative of the peak 
of the spectrum. However, careful inspection of cometary spectra show that lines attributable to 
this band can be seen from approximately 3350 A to 4700 A, although the main part of the band 
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is between 3900 and 4140A with a maximum at 4050A and an additional (but not well identified) 
peak at 4300A. 

The structure of this molecular band is quite complex since the bending frequency of the 
lower state is small and there is a large Renner- Teller splitting. The band has been identified as 
a A 1 ]!^ - X 1 !!^" electronic transition (Gausset et al., 1965). R-branch bandheads at 4072A and 
4260A can be assigned to the (1,0,0)-(1,0,0) and (0,0,0)-(l,0,0) bands respectively (Merer, 1967). 
The density of the lines results in a pseudo continuum from C3. 

C3 is a linear, symmetric molecule containing nuclei which are identical in mass. For all £ 
states of the molecule, every second line is missing. In other states, one member of the e//-parity 
doublet of each rotational level is missing, although all rotational levels are represented (Tokaryk 
and Chomiak, 1997). 

Most of our line identifications come from Gausset et al. (1965), with additional identifica- 
tions from Merer (1967) and Tokaryk and Chomiak (1997). Tokaryk and Chomiak argue convinc- 
ingly that the (0,2,0)^(0,0,0) and (0,2,0)-(0,2,0) rotational assignments of Gausset et al. are in 
error. Thus, in the wavelength region from 3990-4019A we have preferentially used the Tokaryk 
and Chomiak assignments. However, Tokaryk and Chomiak do not report rotational assignments 
for these bands outside this wavelength region, so we have continued to use the Gausset et al. 
assignments for these two bands at other wavelengths. 

The combination of the complexity of the spectrum and the difficulty of studying this molecule 
in the laboratory has meant that the line identifications are far from complete and may not always 
be accurate. This is readily apparent at 3949A where we observe an obvious C3 bandhead which is 
not in the line lists cited here. The lines which are identified in this region arise from the (0,4,0) - 
(0,2,0) band Q-branch, but the line identifications of Gausset et al. do not include any assignments 
for the R- or P-branches. Balfour et al. (1994) have studied C3 at these wavelengths and do indeed 
identify this bandhead as the R-branch of the (0,4,0) - (0,2,0) band and identify other lines in the 
region as from the P-branch. However, the Q-branch assignments of Balfour et al. disagree with 
those of Gausset et al. Tokaryk and Chomiak have pointed out that the Balfour et al. spectra 
are at much lower resolution than the other spectra and that Balfour et al. may not have been 
able to disentangle the superimposed bands. We have no basis for choosing between the various 
assignments. We have chosen to neglect the Balfour et al. assignments entirely based on the 
arguments of Tokaryk and Chomiak. However, it is reasonably certain that the undentified lines 
near 3949A are from the R-band of the (0,4,0) - (0,2,0) band. 

We do not resolve cleanly the C3 lines in our spectra. The outflow of the gas in the coma 
broadens the line sufficiently that even at 3x the resolution with R= 180,000, a resolving power we 
used with comets Hyakutake and Hale-Bopp, the C3 is not resolved in the densest regions. Thus, 
we have not always been able to definitively match a line list with our spectra. In the densest 
regions, we have just assumed that all known lines were observed. There are some known lines (e.g. 
in the (0,0,0)-(0,4,0) or (0,0,0)-0,6,0) bands) which we definitely do not observe in our spectra. 
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3.3 CH 

A wonderful new, inclusive line list has recently become available for studies of CH. This is the 
SCAN-CH list of J0rgensen et al. (1996). This list includes transitions of ^C 1 !! and ^C 1 !! for the 
infrared transition X 2 II - X 2 II and three electronic transitions A 2 A - X 2 n, B 2 S _ - X 2 II, and 
C 2 S + - X 2 n. A total of 112,821 lines are given for each isotope! We identified low-J lines of the 
0-0 band of the B-X transition around 3900A (some of which can be confused with isotopic lines 
of CN) and low-J lines of the 0-0 and 1-1 bands of the A-X transition around 4300 A. This new 
SCAN-CH line list includes the satellite transitions of the CH bands, which are observed in our 
spectra. 

3.4 CN 

Great effort has been expended to understand the CN molecule because it is a significant source 
of opacity in some stars. Thus, there are CN line lists which are much more inclusive than what is 
needed for cometary studies, where depopulation of the upper J-levels is relatively rapid. 

Two different electronic band systems are seen in cometary spectra, the so-called violet 
system, or the B 2 S + - X 2 S + band, and the so-called red system, or the A 2 n - X 2 S + band. 

For the violet system, we used the line list of Kurucz (1995, CDRom 18) which contains 
about 350,000 lines. Lines of both 12 C 14 N and 13 C N are included with J up to 150.5. For our 
purposes, we have only used the 12 C 14 N wavelengths. Indeed, some of our "unidentified" lines 
may be lines of 13 C 14 N. Since this band is a S-S transition, there are P- and R-branches but no 
Q-branch. 

For the red system, we used the line list of Davis and Phillips (1963). J0rgensen and Larsson 
(1990) recently calculated CN red line lists with up to 1 million lines for 12 C 14 N. However, the 
J0rgensen amd Larsson work is most concerned with completeness for line blanketing and is not 
derived from measured line positions. Thus, some of the line center uncertainties in this list are 
higher than we preferred, which is why we used the older reference. 

3.5 NH 2 

Throughout much of the red region of a cometary spectrum, there are scattered emission lines 
which are attributable to NH2. The transitions responsible for these lines are the A 2 Ai - X 2 Bi 
electronic transition along with transitions between the high vibronic levels and the ground state 
of the X 2 Bi electronic band (e.g. (0,13,0)X 2 B x - (0,0,0)X 2 Bi). The spectrum of the NH 2 band is 
quite complex and irregular because the upper level is linear and the lower level is bent at an angle 
of 103°. Also, perturbations are important, as is the Renner- Teller effect. 

The traditional reference for cometary NH 2 line lists is Dressier and Ramsay (1959). While 
this reference is moderately complete, it has been superseded in subsequent years by several other 
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laboratory analyses. Most important in newer studies is the inclusion of more complete perturbation 
models and higher quality laboratory spectra. For our line lists, we have used the data of Ross et al. 
(1988) for the wavelength range 5400 - 6800A. Redward of 7000A we used the data of Johns et al. 
(1976) and of Huet et al. (1996), with updates by Hadj Bachir et al. (1999). From 6800-7000A (the 
(0,2,0)-(0,0,0) and (0,l,0)-(0,0,0) bands), the data of Dressier and Ramsay and Johns et al. (1976) 
are most complete, so we use these two lists preferentially in that region. For regions blueward of 
5400A, the Dressier and Ramsay wavelengths are really the only data in existence; we used their 
wavelengths in this blue region. 

Since NH2 has the perturbations and the bent state, modern works prefer to use the bent 
notation for the vibronic structure. This is in contrast with the work of Dressier and Ramsay and 
Johns et al, in which the vibronic states were denoted in the linear notation. In our line lists, 
we have converted all of the vibronic notation from the linear to the bent notation. In the linear 
notation, sub-bands were denoted by E, II, A, <E>, and T (for K= 0, ±1, ±2, ±3, ±4 respectively). 
To convert from a linear to a bent notation, one uses the formula 

v lin = 2v bent + \l\ 

where I = K ± A. Thus, the old notation for the band around 6500A would have been (0,8,0) II or 
(0,8,0) <I>, while now those bands are denoted (0,3,0) and (0,2,0) respectively [in these examples the 
lower level is assumed to be (0,0,0)]. Despite the fact that the linear notation has been exclusively 
in use by the comet community, we feel it is more appropriate to adopt the same notation as the 
physicists who measure these molecules. 

3.6 CH+ 

The CH + lines which are seen in cometary spectra are from a A 1 !! — X 1 E + transition which 
contains P, Q and R branches. The line list used for these bands was from Douglas and Herzberg 
(1942). The bands degrade towards the red. We only see lines of the 0-0 branch with relatively low 
J-values. While we see very few CH + lines in our spectra, they are visible in both the optocenter 
and the tail spectra. 

3.7 CO+ 

There are three band systems which can be attributed to the CO + molecule. These are the first 
negative carbon bands (B E - X 2 E + ), which occur in the UV, the comet-tail bands (A 2 II - X 2 E+), 
occurring in the violet through red, and the Baldet-Johnson bands (B 2 E + - A 2 II) also in the violet. 
Observations in the optical bandpass do not include the first negative carbon bands (however, UV 
observers can find wavelengths for these lines in Rao (1950a). Rao (1950b) showed that previous 
optical band assignments were in error and that the Baldet-Johnson system was much weaker than 
the comet-tail bands. Thus, we adopted the band assignments of Rao (1950b) for the comet-tail 
bands and neglected the Baldet-Johnson system entirely. 
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For our comet-tail line list, we used four references: Coster et al. (1932), Schmid and Gero 
(1933), Rao (1950b) and Haridass et al. (2000). We updated the band designations of the first 
two references with the designations in Table 2 of Rao and also corrected the values for N in these 
papers as outlined by Rao. Unlike for CH + and R^O" 1- , the CO + lines could not be detected in the 
optocenter spectra. These were entirely found in the tail spectrum. Indeed, we detected only the 
2-0 band of CO + in the de Vico tail spectrum. 

It is unusual to detect the 2-0 band without the 3-0 band since the 3-0 band is normally the 
strongest CO + band in cometary spectra. The ion tail of de Vico was generally quite narrow so it 
is possible that our "tail" spectrum missed the strongest emissions, accounting for the weakness of 
the CO + emissions. 

3.8 H 2 0+ 

With H2O being the dominant ice in the cometary nucleus, it is not surprising to see H20 + emission 
in the ion tail spectrum of a comet. The electronic transition is A 2 Ai - X 2 Bi, which is observed 
from 4000-7500A. Much of the structure of the H20 + molecule is similar to that of the NH2 
molecule, including the strong Renner- Teller effect. 

The most complete reference for the wavelengths of the H20 + lines is Lew (1976). In this 
paper, the bands are denoted in the linear notation (see the section on NH2) whereas the bent 
notation is considered correct by current standards. Therefore, we have converted the H20 + bands 
to the bent notation in a manner which is identical to the NH2. 

As with the CH + , the H20 + lines are apparent in both the optocenter and the tail spectra. 
We required that a line be visible in the tail spectrum before we attributed it to H20 + . 

3.9 Night Sky Emission Lines 

Finally, even at the high spectral resolving power of our spectra, coupled with the small entrance 
slit, one does end up observing emission lines produced in the Earth's atmosphere. These lines 
must be accounted for when identifying lines. Fortunately, a high-resolution atlas of night sky 
emission lines has recently been compiled by Osterbrock et al. (1996) for just this purpose. The 
dominant molecular emission comes from OH, with contributions of O2 and some atomic lines. We 
have applied a Doppler shift equal to —A to the Osterbrock et al. wavelengths to compensate 
for correcting the cometary spectra to the laboratory rest frame. We confirmed that the correct 
Doppler shift was applied by checking the strong night sky atomic lines (Na D, O ( X D) and H a). 

In the lists of identified lines which follow, we do not include the night sky lines, though they 
are represented on our plots. We have drawn no conclusions about their existence in our spectra, 
assuming they are all present. Instead, an otherwise unidentified line which is coincident with a 
night sky line is considered to be that night sky line without taking into account the line's strength. 
Some unidentified cometary and night sky lines could, however, be coincident. 



9 



3.10 Additional Molecules 

The molecules which have been identified above are the same ones which Brown et al. (1996) found 
in their spectra and which previous authors have identified. Brown et al. also investigated a range 
of other molecules which did not show up in their spectra and we have assumed they are not in our 
spectra either. 

There are some potential molecules which we checked but our investigations proved incon- 
clusive. One is H2O, which is the dominant molecule in cometary nuclei and which has many 
transitions in our spectral region. The wavelengths of these lines are well known. However, though 
we started out examining our spectra for H2O lines, we quickly realized that it would be implausible 
that we could detect these lines in our spectra when the Earth's atmosphere is so full of H2O. 

Another potential molecule in cometary spectra is N2. Again, the atmosphere of the Earth 
makes this molecule impossible to study, so we did not search for it. 

One molecule for which we searched was . Cochran et al. (2000) showed that this molecule 
was not present in spectra of de Vico or Hale-Bopp and placed stringent upper limits on its presence. 

4 Line Identifications 

Once the molecular line lists were compiled, it was necessary to examine all of the spectral orders 
and determine which lines we detected. First, we scrutinized our spectra for wavelength coincidences 
of the emission features and lines in the line lists. The strongest lines are quite easy to see in the 
spectra, but the weakest lines needed verification. We utilized spectra from both nights to test the 
veracity of a detection of a weak line. We required that the line be in both night's spectra before 
we would categorize the feature as a detection. Our very high signal/noise, coupled with having 
multiple spectra, allowed for the detection of some very weak features. 

With the density of the lines in some of these bands and in our spectra, it is expected that 
there will be some coincidences which are accidental. When a line from the line list was coincident 
with an emission line, we checked to determine whether detection of the line made sense. This 
criterion took the form of requiring that the lines originating from lower energy levels be observed 
before believing that lines from higher energy transitions were detected. This is equivalent to 
assuming that the upper level populations can be described by some sort of effective Boltzmann 
temperature, which is not necessarily a physical temperature. Indeed, because C2 is a homonuclear 
molecule, with no allowed dipole transition, high-J levels can be quite strong because these levels 
get "pumped up". This criterion was applied on a band-by-band basis. Thus, there were times 
when an "obvious" detection was deemed accidental because lines with lower excitation energies 
for the same band were not seen. We would then eliminate these accidental detections from our 
list. 

Conversely, there are some times when we included a line, though it is not very obvious, 
because it is more favorable energetically than lines which were detected. A possible cause for the 
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absence of an expected line is coincidence with a solar absorption feature. 

For the ions, we checked the tail spectrum as well as the optocenter spectra in order to 
determine if a line was ionic in nature. A simple coincidence of the ionic species wavelength to 
a feature in the optocenter spectra was not considered sufficient evidence that a line was an ion. 
The line also had to be present in the tail spectrum. As noted above, the CH + and H20 + lines 
were observed in both the optocenter and tail spectra, but the CO + was only observed in the 
tail spectrum. It was considered sufficient to identify an ionic line if it occurred solely in the tail 
spectrum. 

We compared our tail spectrum of de Vico with the list of identified and unidentified species 
in the spectrum of the tail of comet C/1996B2 (Hyakutake) (Wyckoff et al. 1999). Some of 
Wyckoff et al. 's H^O" 1 " lines are seen in our spectra, while others are not, and we detect many 
H20 + lines that they did not see (however, our data are at much higher spectral resolution and 
so exact matches in wavelength with their identifications are difficult). We cannot confirm most 
of their lines; though a few of their lines other than H20 + might be coincident with features in 
the de Vico tail spectrum, these features are generally also coincident with neutral lines identified 
in the optocenter spectra. The Hyakutake tail spectrum of Wyckoff et al. was obtained at larger 
distances from the nucleus than was our de Vico tail spectrum, which may explain why we do not 
see the same lines. In addition, with the high airmass of observations and a relatively small slit on 
our instrument, we could have missed the densest part of the ion tail (we do see some ions though). 
However, the explanation could as easily be that de Vico does not have a very strong tail spectrum. 

In the blue part of the spectrum, there was an increasing wavelength overlap between the 
echelle orders as the wavelength decreased (higher order number). This allowed for the comparison 
of many lines in two different echelle orders. In general, we found the agreement excellent, with 
wavelength matches of ~ 0.01 A or better. However, for the bluest orders, the spectral overlap 
regions were falling well off the blaze of the grating and, consequently, became noisy in the bluest 
quarter of these orders. Therefore, we gave priority for identification of lines in the overlap regions 
to the red end of one order over the blue end of the next lowest (and therefore redder) numbered 
order. 

One interesting result of our identifications comes from the C2 Swan bands, particularly the 
Av = —1, 0, andl bands. We observe lines with very high J-values, up to J=109! However, the 
lowest few J-value lines are missing or very weak. This is probably indicative of two different 
rotational temperatures for C2, though we have not modeled these lines yet to demonstrate this. 
Such an effect was also seen by Lambert et al. (1990) for comet Halley. 

For the CN red system, we found that for a particular J level, the P12 and R12 lines were 
generally quite weak when compared with the P, R or Q lines of the same J level. Inspection of 
the line strengths in the laboratory work of Davis and Phillips (1963) shows that this is generally 
seen for this molecule and is normal. 

In addition to the molecular emission lines, we found three atomic emission lines in the 
spectrum. These are the red metastable oxygen doublet, O (*D) 6300A and 6364A, and H a. The 
green oxygen line at 5577A is cloaked by the 1-2 C2 emission band; the Na "D" doublet is in the 
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wavelength region of the C2 Af = —2 bands. However, at most, there is just a weak hint of the 
cometary Na features at appropriate wavelengths. The weakness of any Na emission in de Vico is 
quite unusual in a comet at de Vico's heliocentric distance. 

Finally, when we had exhausted our line lists, we still had a great many emission lines in our 
spectra with no identifications. We have flagged these in our line lists with the tag "Unid" next 
to the measured wavelengths and computed vacuum frequencies. A few of these "unidentified" 
lines are probably isotopic lines, e.g. 13 CN vs. 12 CN. Isotopic lines are quite weak but we do 
know of several in our spectra. However, we have chosen to ignore them for now, because proper 
identification takes modeling the band in question, which is beyond the scope of this atlas. 

Figure [l| shows a plot of the blue half of order 49, covering the wavelength range from 6957- 
7014A. In this fi gure, we show spectra from both nights, offset from one another for clarity. The 
spectra are plotted with an expanded y axis in order to illustrate the weaker features. The spectra 
are plotted in three panels to allow for inspection of the features. Plotted below the spectra are 
tick marks for all of the identified lines, in this case lines of CN, C2, NH2 and H20 + . Plotted above 
the spectra are tick marks for the night sky OH (Doppler shifted to the appropriate rest frame for 
the comet) and the position of lines we could not identify. Dashed vertical lines mark the position 
of the features in the spectra. Examination of Fig. |l] shows that there are a substantial number of 
quite strong lines which are unidentified, including lines with more than 20,000 counts! 

In contrast, Figure |2| shows order 77 from 4 October. There are obvious emission features 
(and solar absorptions) in this spectrum. Indeed, we have located 64 lines in this spectral region 
which are visible in spectra from both nights. However, we can not identify the molecule for any 
of these 64 features. This order is an excellent example of the desperate need for more laboratory 
studies. 

Figure || shows a plot of all of our detected lines, by molecule, as a function of wavelength. 
Unidentified features are also marked. The gaps which can be seen in the red for molecules such 
as C2 and CN are the result of the interorder gaps in our spectra. Molecules such as C2 and NH2 
can be found at almost any optical wavelength, while C3 and CH are confined to the blue. 



Tabic III gives a representative listing of the identifications for the wavelength region from 



6957-6976A, the spectral region shown in the top panel of Fig. [I]. A table similar to Table [II for 
each of the molecules shown in Figure ^ is available in digital (ASCII text) formP]. These tables 
contain the wavelength and vacuum frequency of the lines, along with the band designations and 
transition. A table of the unidentified lines is also available in digital form. We also will make 
available a table of all of the identified and unidentified features, sorted by wavelength. 



As with Table [II, the digital tables will be entirely in ASCII characters and will not use any 
sub- or super-scripts; they can be utilized on any computer without having to deal with embedded 
special characters. The tables are in fixed format. For most molecules, the electronic transitions 
are listed by their letter designations (A, B, X) instead of their term designations. For C2, we use 
the more common names of Swan and Phillips, where the Swan band is the d 3 n„ - a 3 n u band 



1 These tables and the plots mentioned below may be obtained from the web site (currently pdssbn.astro.umd.edu) 
of the Small Bodies Node of the NASA Planetary Data System 
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and the Phillips is the A H u - X £jj" band. As vibrational transitions carry designations with 
different forms for diatomic and polyatomic molecules, those transitions are listed in the manner 
appropriate for each molecule. In all cases, the upper state is listed first, then the lower state. For 
the rotational transitions, the designations are very dependent on the individual molecule, since 
some just have P, Q and R branches, some, such as CN A-X, have branches such as P12, and some, 
such as CH, have branches such as ^Q2i/e (which is listed as " A RQ21fe" in plain ASCII). Molecules 
such as NH2 and H20 + have rotational states which are designated with three quantum numbers 
of the form N' K , K , — N'J^,,^,,. The NH2 and H20 + rotational states are separated into Fi and F2 
components. We have not included those designations in our tables. 

Sometimes more than one transition occurs at the same frequency. We list all rotational 
transitions from the same vibrational and electronic band on the same line in the tables. When 
there are concurrent transitions from different vibrational bands, we list the two different band 
transitions on separate lines of the table. 

We chose to list only the wavelength and frequency of the unidentified lines and not their 
strengths. There were several reasons behind this decision. First, we have no absolute flux calibra- 
tion of these coude spectra so there is no way to calibrate the blaze function. Additionally, we have 
not attempted to remove the atmospheric extinction and the comet was observed at an airmass 
of 2-3 airmasses. This means that relative strengths within an order are probably accurate, but 
order-to-order comparisons are less meaningful. We have made no attempt to model the Swings 
effect for any molecule. Without such a model, any line could be a different strength than expected 
from laboratory measurements. Also, the lines are often highly blended. Thus, to determine line 
strengths a full deconvolution model must be computed. For the over 4000 unidentified lines this 
is a huge undertaking. Such an exercise should be done for individual studies, but are not realistic 
for this atlas. 

In addition to the tables of identified and unidentified features, we are making available 
PostScript files of plots of our spectra in two formats. The first format is plots similar to Fig. [l|. 
For each order, there are two plots, each with 3 panels, showing two spectra and the identified and 
unidentified lines. These are all plotted with expanded y scales to show the weaker features. The 
features are marked in color, as with Fig. [l]. The colors of the markers remain the same for each 
molecule across all plots. That is, there is space and a color allocated for each molecule even when 
that molecule is not present in a particular order. Included in these plots are the ionic species CH + 
and H20 + , since they are visible in the optocenter spectra. The CO + only shows up in order 81 
in the tail spectrum. This spectrum is shown in Fig. 2 of Cochran et al. (2000) with the CO + and 
CH+ identified. 

The second form of plot which we are making available is the spectra without expansion of 
the y axis. These plots take the form of four panel plots with two orders per plot (each of two 
panels). Each panel is scaled in y to fit the maximum strength line in that panel. These plots are 
a valuable adjunct to the identification plots because they show the relative strength of the strong 
and weak lines. They also allow for the inspection of the telluric absorption features which plague 
some orders. 

It is quite likely that there are a few real features which we have missed in our tabulations 



13 



because the sheer quantity of data makes such misses inevitable. Indeed, every time we have 
inspected the plots of identifications, we have found an additional line or two. There are some 
lines which we may have located which may be less believable to a reader. Sometimes our choice 
to include these lines in our lists comes from our ability to manipulate the spectra digitally and 
to change the scaling to view a feature better; sometimes we may have overactive imaginations! 
There are times when an identified molecular line seems slightly off in wavelength from the peak of 
a feature. We have assumed that it is possible that the wavelengths of the laboratory line lists are 
not always entirely accurate. Additionally, in regions where lines are blended, or in the regions of 
strong telluric absorption features, the observed line center may be shifted from a laboratory value. 
This is true also for our measured positions of unidentified features. Therefore, it is imperative 
that a user of these data inspect the full scale plots for context and for absorptions and that they 
examine the feature plots for blending. Users who wish to model the unidentified features should 
contact us for access to the data, so that the lines can be properly deblended. 

5 Summary 

In 1995, we obtained high spectral resolution observations of the dust-poor comet 122P/de Vico 
at McDonald Observatory. We have compiled many laboratory studies of molecules thought to be 
present in the spectra of comets and have used these to identify 12,219 lines in the spectrum of 
de Vico. In addition, we found 4,055 lines we could not attribute to a particular molecule. We 
found emission lines of five neutral radicals: 6862 were of C2, 1282 of C3, 1167 of CN, 169 of CH, 
and 2569 of NH2. We also identified three atomic lines. We observed three ionic species: 129 lines 
of H20 + , 11 lines of CH + and 27 lines of CO + . The CO + can only be seen in the spectrum obtained 
100 arcsec tailward, while the H2O 4 " and CH + can be seen in both the tail and optocenter spectra. 

Our lists of identified and unidentified emission lines are available in digital format from the 
Planetary Data System Small Bodies Node. In addition, we are making available annotated plots 
of high resolution spectra of a comet, covering almost the complete optical bandpass. The line lists, 
coupled with the plots, should prove to be an invaluable tool for future studies of high resolution 
cometary spectra. Only by first identifying the molecules which are present in a spectrum can the 
important chemical analyses be performed. 
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Table I: Log of Observations 
Start 





Date 


Time 


Exposure 


A 


A 


Rh 


Rh 


File # 


(UT) 


(UT) 


(sec) 


(AU) 


(km/sec) 


(AU) 


(km/sec) 


RV23200 


03 Oct 95 


11:37:30 


600 


1.00 


-14.3 


0.66 


-2.9 


RV23201 


03 Oct 95 


11:57:17 


600 


1.00 


-14.3 


0.66 


-2.9 


RV23257 


04 Oct 95 


11:10:49 


1500 


0.99 


-12.9 


0.66 


-1.7 


RV23258" 


04 Oct 95 


11:42:21 


1200 


0.99 


-12.9 


0.66 


-1.7 


Note: 


a) 100 arcsec 


tail ward 













Table II: Wavelengths Observed 



Order 


Start 


End 


Order 


Start 


End 


Order 


Start 


End 




(A) 


(A) 




(A) 


(A) 




(A) 


(A) 


89 


3829.77 


3893.62 


70 


4869.76 


4950.90 


51 


6684.03 


6795.33 


88 


3873.32 


3937.89 


69 


4940.35 


5022.66 


50 


6817.70 


6931.23 


87 


3917.88 


3983.19 


68 


5013.01 


5096.53 


49 


6956.83 


7072.67 


86 


3963.46 


4029.53 


67 


5087.84 


5172.60 


48 


7101.75 


7220.00 


85 


4010.12 


4076.96 


66 


5164.94 


5250.98 


47 


7252.84 


7373.60 


84 


4057.89 


4125.53 


65 


5244.41 


5331.77 


46 


7410.50 


7533.88 


83 


4106.81 


4175.26 


64 


5326.36 


5415.08 


45 


7575.16 


7701.29 


82 


4156.91 


4226.20 


63 


5410.91 


5501.04 


44 


7747.31 


7876.30 


81 


4208.26 


4278.40 


62 


5498.18 


5589.76 


43 


7927.46 


8059.45 


80 


4260.88 


4331.90 


61 


5588.32 


5681.40 


42 


8116.19 


8251.32 


79 


4314.84 


4386.75 


60 


5681.46 


5776.09 


41 


8314.13 


8452.55 


78 


4370.18 


4443.01 


59 


5777.75 


5873.98 


40 


8521.97 


8663.85 


77 


4426.96 


4500.73 


58 


5877.37 


5975.26 


39 


8740.46 


8885.98 


76 


4485.23 


4559.97 


57 


5980.48 


6080.08 


38 


8970.45 


9119.79 


75 


4545.05 


4620.78 


56 


6087.27 


6188.65 


37 


9212.88 


9366.25 


74 


4606.48 


4683.24 


55 


6197.94 


6301.16 


36 


9468.77 


9626.40 


73 


4669.60 


4747.41 


54 


6312.71 


6417.84 


35 


9739.28 


9901.42 


72 


4734.47 


4813.35 


53 


6431.82 


6538.93 


34 


10025.71 


10192.61 


71 


4801.16 


4881.16 


52 


6555.50 


6664.66 
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Table HI: The Lines of Order 49, panel 1 



V V Oj V OlCIlg Ull 


n rpn 1 1 (*y\ r*\i 
1. 1 c-v^u.C'ii^^y 






Vibrational 2 


T\ of" atinn a 1 


(A) 




A/Tolooi 1 1 

IV-HJlCv/ Lilt 




Transition 


TV fi n q l i~ l on t; 

11 clllolLlUllo 




14d7(J.o97 


H2U+ 


A— A 


(0,2,0) - 


- (0,0,0) 


O 1 o o o o 

3 12-202 


r-r\ r T OA1 

6957.301 


1 A •) t ' <\ A OO 

14369.430 


/"I AT 


A- A 


3- 


-0 


Pl( 2) 


r*c\ r t c* 1 o 

6957.612 


14368.780 




A "V 

A-A 


3- 





no/ o\ 

P2( 9) 


f'C\ r !"7 1701 

6957.731 


14368.540 




A V 

A-A 


3- 


-0 


/ \ -1/ o\ n-i /-|n\ 

C;l( 8; Kl(19) 


6958.122 


1 ^ or"7 Ton 

14367.730 


ATTTO 

IN 112 


A V 

A-A 


(0, 2,0) 


- (0,0,0) 


r 1 r roc 

5 15-505 


6958.228 


1 /i no 

14367.513 


H2U+ 


A V 

A-A 


(0,2,0) - 


- (0,0,0) 


O 1 o o o o 

6 1 2 - 2 (J 2 


6958.640 


1 A^c c z' 1 r r 

14366.655 


T T "J 

Unid 










6958.864 


1 A oca onn 

14366.200 


C2 


ID T, • 1 

rhil 


4- 





Q(18) 


/^oro oo'v 

6958.907 


14366.110 


CM 


A V 

A-A 


3 


-0 


Pl( 3) 


6959.100 


14365.710 


CJN 


A V 

A— A 


3 


-0 


Cj2(14) 


6959.394 


14365.100 


CJN 


A V 

A-A 


3 





Ql( 9) 


6960.54/ 


1 /i 0^0 79n 
14362. M0 


CJN 


A Y 

A— A 


3- 


-0 


\xY\2\)) 


enen one; 

6960.805 


1 /i oeo 1 nn 

14362.190 


CJN 


A V 

A— A 


3- 


-0 


TD 1 / A \ 

Pl( 4) 


6960.850 


1 a oci 1 nn 

14362.100 


ATTJO 
IN 112 


A V 

A— A 


(0, 2,0) 


- (0,0,0) 


C 1 c c o c 

1 5 - U 


6961.318 


1 /i o£; 1 1 on 
14361.130 


CJN 


A Y 

A— A 


3- 


-0 


oi /-i o\ 


6961.720 


1 /i oen 010 

14360.312 


Umd 










ncia~\ nno 

6961.903 


1 /i ocn non 

14359.930 


CJN 


A V 

A-A 


3 


-0 


r , 2{W) 


6962.377 


1 /loco ncn 

14358.950 


C2 


rhil 


4- 





P(14) 


6962.795 


1 /loco non 

14358.090 


CJN 


A V 

A— A 


3 


-0 


ni o/ o\ 00/1 c; 1 ! 

rl2{ 6) yz(lo) 


n ore 

6962.955 


14357.7bU 


/"I AT 

CJN 


A V 

A-A 


3- 


-0 


Til ( C\ 

Pl( 5) 


6963.331 


1435b. 


H2U+ 


A V 

A-A 


(0,2,0) - 


- (0,0,0) 


O 1 1 1 O 1 

2 1 1 — 1 U 1 


6963.496 


1 /lore £!/in 


CJN 


A Y 

A— A 


3- 


-0 


oi ( -\-\\ 

yi(ii) 


6963.569 


1 /l O ET £? /I OO 

14356.490 


/"I AT 

CJN 


A V 

A-A 


3- 


-0 


I ~> 1 /r\i \ 

Wa\2\) 


6963.594 


1 /i occ /in 

14356.440 


ATTJO 
IN 112 


A V 

A-A 


(0, 2,0) 


- (0,0,0) 


A 1 A A O A 

4 1 4 - 4 (J 4 


6963.880 


14355.862 


T T~ ; A 

Umd 










r> r\r' A O OO 

6964.080 


1/iocrr - a a 'i 

14355.447 


TT 'J 

Umd 










f*C\f A OAO 

6964.303 


1 AOtZ A HOA 

14354.980 


ATTJO 
IN 112 


A V 

A-A 


(0, 2,0) 


- (0,0,0) 


111 1 O 1 

1 1 1-1 U 1 


r> r\r' a r no 

6964.590 


1 A O Cf ,4 OO /I 

14354.394 


TT 'J 

Umd 










ccic A onn 

6964.800 


1 /i r nc 1 

14353.961 


tt~ ; J 
Umd 










6965.392 


1 /I O C O 'I A C\ 

14352.740 


/"I AT 

CJN 


A V 

A-A 


3- 


-0 


Pl( 6) 


6965.470 


14352.575 


H2U+ 


A V 

A-A 


(0,2,0) - 


- (0,0,0) 


o 1 1 1 n 1 

2 11-10 1 


Dyoo.yui 


14lo01.DyU 


PINT 


A Y 


3- 


-0 


ni /i 9"\ 


6966.064 


14351.350 


NH2 


A-A 


(0, 2,0) 


- (0,0,0) 


5 14-422 


6966.200 


14351.070 


C2 


Phil 


4- 


-0 


R(26) 


6966.212 


14351.050 


CN 


A-A 


3- 


-0 


P12( 4) 


6966.326 


14350.810 


NH2 


A-A 


(0, 2,0) 


- (0,0,0) 


4 14-404 


6966.452 


14350.550 


CN 


A-A 


3- 


-0 


P2(ll) 


6966.637 


14350.170 


NH2 


A-A 


(0, 2,0) 


- (0,0,0) 


2 12-202 


6966.700 


14350.040 


NH2 


A-A 


(0, 2,0) 


- (0,0,0) 


3 13-303 


6966.719 


14350.000 


CN 


A-A 


3- 


-0 


Q2(16) 


6967.060 


14349.305 


Unid 
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Table (cont.) 



Wavelength 


Frequency 




Electronic 1 


Vibrational 2 


Rotational 3 


(A) 


(cm' 1 ) 


Molecule 


Transition 


Transition 


Transitions 


6967.380 


14348.646 


Unid 








6967.730 


14347.909 


Unid 








6967.904 


14347.560 


NH2 


A-X 


(0, 2,0) - (0,0,0) 


5 14-422 


6968.061 


14347.240 


CN 


A-X 


3-0 


Pl(7) 


6968.260 


14346.832 


Unid 








6968.584 


14346.160 


CN 


A-X 


3-0 


Ql(13) 


6968.760 


14345.798 


Unid 








6969.124 


14345.050 


C2 


Phil 


4-0 


Q(20) 


6969.396 


14344.490 


NH2 


A-X 


(0, 2,0) - (0,0,0) 


3 12-220 


6969.993 


14343.260 


NH2 


A-X 


(0, 2,0) - (0,0,0) 


3 13-303 


6970.001 


14343.240 


CN 


A-X 


3-0 


P12( 5) 


6970.152 


14342.934 


H20+ 


A-X 


(0,2,0) - (0,0,0) 


110-000 


6970.499 


14342.220 


NH2 


A-X 


(0, 2,0) - (0,0,0) 


111-10 1 


6970.846 


14341.500 


CN 


A-X 


3-0 


Q2(17) 


6970.892 


14341.410 


NH2 


A-X 


(0, 2,0) - (0,0,0) 


2 12-202 


6971.032 


14341.120 


CN 


A-X 


3-0 


Pl(8) 


6971.181 


14340.820 


CN 


A-X 


3-0 


P2(12) 


6971.513 


14340.130 


CN 


A-X 


3-0 


Ql(14) 


6971.900 
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NOTES 



1: For C 2 , "Phil" is the Phillips band (A - X 

2: The designation of the vibrational transition is different for diatomic and 

polyatomic molecules. 
3: The rotational transition designations depend on the molecule. 

For CN, more than one transition may be listed per line if they occur at 

the same wavelength. 

For NH2, rotational quantum numbers are listed as N' K , K , — N'^„ K „. 
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6 Figure Captions 



Figure H]: A representative spectral order is shown with an expanded y scale in order to demon- 
strate the detection of weak features. Two spectra are plotted. The bottom, black spectrum, 
is RV23257, from 4 October 1995. The top, blue spectrum, is RV23201, from 3 October 1995. 
The blue spectrum is offset vertically for clarity. Marked on the plot are all of the identified and 
unidentified features in the spectrum. 

Figure |2|: Spectral order 77 shows many strong emission features but we can not identify the 
molecule responsible for any of these lines. There are 64 emission lines in this spectral region. 

Figure |||: The wavelengths of the detected lines of the various molecules and atoms we detected 
are plotted. The gaps at the red end for molecules such as C2 are caused by the interorder gaps in 
our spectra. The wavelengths where we have unidentified lines are also marked. 
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122P/deVico Order 49 (part a) 
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Figure 1: Cochran and Cochran 2001 
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122P/de Vico -- Order 77 (RV23257) 
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Figure 2: Cochran and Cochran 2001 
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